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Abstract. Since the Maxwell theory of electromagnetic phenomena is a 
gauge theory, it is quite important to evaluate the zero-point energy of the 
quantized electromagnetic field by a careful assignment of boundary condi- 
tions on the potential and on the ghost fields. Recent work by the authors 
has shown that, for a perfectly conducting spherical shell, it is precisely the 
contribution of longitudinal and normal modes of the potential which enables 
one to reproduce the result first due to Boyer. This is obtained provided that 
one works with the Lorenz gauge-averaging functional, and with the help of 
the Feynman choice for a dimensionless gauge parameter. For arbitrary val- 
ues of the gauge parameter, however, covariant and non-covariant gauges 
lead to an entangled system of three eigenvalue equations. Such a problem 
is crucial both for the foundations and for the applications of quantum field 
theory. 
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Since Casimir produced his remarkable calculation of zero-point energy 
of the electromagnetic field for the case of two perfectly conducting plates 
of area A and separation d: AE = — ^^|f Jr, there has been an increasing 
series of efforts in the literature to evaluate suitable differences in zero-point 
energies for various sorts of geometries and in a variety of media. Since this 
is a research field where the theoretical predictions have been tested against 
observations 0, it is clear why the investigation of Casimir energies has 
attracted efforts over half a century, with even better perspectives for the 
years to come. Theoretical physicists, however, who are also interested in 
the general structures, may approach the Casimir energy calculation with 
the aim of re-deriving a non-trivial prediction in a way which is viewed as 
more fundamental from the point of view of general principles. In particu- 
lar, we refer here to the path-integral quantization of gauge theories. As is 
well known, on using the Faddeev-Popov formalism, one performs Gaussian 
averages over gauge functionals x^i by adding to the original Lagrangian a 
gauge-averaging term x^P^X^ '> where /3 Mi , is any constant invertible matrix. 
This turns the original operator on field perturbations into a new operator 
which has the advantage of being non-degenerate. The result of a one-loop 
calculation is expected to be x~ an d /3-independent, although no rigorous 
proof exists on manifolds with boundary. In particular, for Maxwell theory, 
(3^ reduces to a 1 x 1 matrix, i.e. a real-valued parameter, and x^ reduces 
to the familiar covariant or non-covariant gauges for Maxwell theory, e.g. 
Lorenz, Coulomb, axial. 
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Our paper describes recent work by the authors || on the application 
of the Lorenz and axial gauges to the evaluation of the zero-point energy of 
a perfectly conducting spherical shell. The gauge-invariant boundary con- 
ditions applied by Boyer in his seminal paper [ffl] require that tangential 
components of the electric field should vanish on a two-sphere of radius R. 
More precisely, on denoting by r, 9, ip the spherical coordinates, which are 
appropriate for the analysis of a conducting spherical shell, one writes 

[Ee]aM = 0, (1) 

[E v ] gM = 0. (2) 

The modes of the electromagnetic field are then split into transverse electric 
(TE) or magnetic multipole: f-E = E r = 0, and transverse magnetic (TM) 
or electric multipole: f-B = B r = 0. On using the standard notation for 
spherical Bessel functions, one finds that, in the TE case, Eqs. (0) and (@) 
lead to 

A l j l (kr)+B l n l (kr)] m = 0, (3) 

while in the TM case the vanishing of Eg and E v at the boundary implies 
that 

= 0. (4) 



— (riQjtikr) + Dim(kr))) 



J dM 

Of course, if the background includes the point r = 0, the coefficients B\ and 
Di should be set to zero for all values of I to obtain a regular solution. Such 
a singularity is instead avoided if one studies the annular region in between 
two concentric spheres. 



On the other hand, if one follows a path-integral approach to the quan- 
tization of Maxwell theory, it is well known that the second-order operator 
acting on A^ perturbations when the Lorenz gauge- averaging functional is 
chosen turns out to be, in a flat background, 

p^ = -^D + (i-^)v m v,. (5) 

With a standard notation, g is the background metric, a is a dimensionless 
parameter, []] is the D'Alembert operator 

d 2 

where A is the Laplace operator (in our problem, A is considered on a disk). 
This split of the Q operator leads eventually to the eigenvalue equations for 
the Laplace operator acting on the temporal, normal and tangential compo- 
nents of the potential (see below). At this stage, the potential, with its gauge 
transformations 

£ A^ = A^ + V M £, (6) 

is viewed as a more fundamental object. Some care, however, is then neces- 
sary to ensure gauge invariance of the whole set of boundary conditions. For 
example, if one imposes the boundary conditions 

[At] dM = 0, (7) 

[A e ]dM = 0, (8) 
[A^qm = 0, (9) 
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this is enough to ensure that the boundary conditions flip and (0) hold. 
However, within this framework, one has to impose yet another condition. 
By virtue of Eq. (|j), the desired boundary condition involves the gauge 
function: 

[e} 9 M = 0. (10) 

Equation fllOD ensures that the boundary conditions are preserved 

under the gauge transformations (@) . What happens is that e is expanded in 
harmonics on the two-sphere according to the relation 



oo I 

e(t,r,9,<p) = J2 E ei(r)*im(Me 

1=0 m=-l 

After a gauge transformation, one finds from (H) 



iut 



(11) 





[M 


8M — \St e \dM-, 


(12) 




[As 


8M = [V 'o£\dM, 


(13) 


[ £ ^4 v ]om — 




8M = [V^ejaM, 


(14) 


one has 










8M z 


= iiv[e] dM , 


(15) 



oo I 

Yl Yl £i(j-)Yi mfi {6,ip)e 

1=0 m=-l 



LLUt 



OO I 

E E £ l( r ) Y lmA e > f) e 
1=0 m=-l 



iujt 



dM 



dM 



(16) 
(17) 



Thus, if £i(r) is set to zero at the boundary V/, the right-hand sides of fll5|)- 
(|T3) vanish at dM. But 

[e { (r)W = ovz 



is precisely the condition which ensures the vanishing of e at DM, and the 
proof of our statement is completed. 

At this stage, the only boundary condition on A r whose preservation 
under the transformation (Q) is again guaranteed by Eq. ([R]) is the vanishing 
of the gauge- averaging functional <&(A) at the boundary: 

[$(A)] BM = 0. (18) 

On choosing the Lorenz term 3>l(A) = VM^, Eq. (|I8"D leads to 

= 0. (19) 

dM 

It is only upon considering the joint effect of Eqs. flTD-QlOD, (|18[) and flT9|) that 
the whole set of boundary conditions becomes gauge-invariant. This scheme 
is also BRST- invariant. In the following calculations, it will be enough to 
consider a real-valued gauge function obeying Eq. fllOD at the boundary, 
and then multiply the resulting contribution to the zero-point energy by —2, 
bearing in mind the fermionic nature of ghost fields for the electromagnetic 
field. 

In our problem, the electromagnetic potential has a temporal component 
A t , a normal component A r , and tangential components Ak (hereafter, k 
refers to 9 and ip) . They can all be expanded in harmonics on the two-sphere, 
according to the standard formulae 

oo I 

A t (t, r,9,<f)=J2 E <r)Yi m (0, v)e luJt , (20) 

1=0 m=-l 
oo I 

A r (t, r,0,^) = E E bi{r)Yi m {0, tp)e iut , (21) 



9A r 2 A 
or r 



oo / 



(22) 



A k (t, r,6,cp) = J2H [ci{r)d k Y lm {6, <p) + T^e^Y^e, <p\ 

1=1 m=—l 

This means that we are performing a Fourier analysis of the components 
of the electromagnetic potential. The two terms in square brackets of Eq. 
( |22|) refer to longitudinal and transverse modes, respectively. On setting 
a = 1 in Eq. (||), the evaluation of basis functions for electromagnetic 
perturbations can be performed after studying the action of the operator 
—g^ v \Z\ on the components (pO|)-(p^). For this purpose, we perform the 
analytic continuation w — > iM, which makes it possible to express the basis 
functions in terms of modified Bessel functions (of course, one could work 
equally well with u, which leads instead to ordinary Bessel functions). For 
the temporal component, one deals with the Laplacian acting on a scalar 
field on the three-dimensional disk: 



(( 3 ) A A) t 



d 2 A t 



2d A t 



dr 2 r dr 
which leads to the eigenvalue equation 



d 2 + 2d_ 1(1 + 1) 



4 (2) AA) ( 



(23) 



ai = M 2 a\. 



(24) 



dr 2 r dr 

Hereafter, we omit for simplicity any subscript for M 2 , and the notation for 
the modes will make it sufficiently clear which spectrum is studied. The 
solution of Eq. (p4|) which is regular at r = is thus found to be 



ai (r) = —I l+1/2 (Mr) 



(25) 



up to an unessential multiplicative constant. 



The action of on the component A r normal to the two-sphere is 



; {3) a A) r 



where the stroke | denotes two-dimensional covariant differentiation on a two- 
sphere of unit radius. Last, the Laplacian on tangential components takes 
the form 

d 2 A, 



(® A A)i 



dr 2 



+ \( {2) AA) k -\A k + -d, 



T kA r 



(27) 



By virtue of the expansions ( |2"T| ) and 
erties of spherical harmonics, Eqs. (| 
equation 

" d 2 1(1 + 1) 
dr 2 r 2 



, jointly with the standard prop- 
and (|2~7D lead to the eigenvalue 



(28) 



Ti = M 2 Ti 



for transverse modes, jointly with entangled eigenvalue equations for normal 
and longitudinal modes (here I > 1): 

(1(1 + 1) + 2) 



d 2 2d 
dr 2 r dr 



dr 2 



1(1 + 1) 



21(1 + 1) 
h H 5 c, 



c, + -6, = M c;. 
r 



(29) 



(30) 



The mode &o( r ) is instead decoupled, and is proportional to /3/2 J^ fr - ) in the 
interior problem. It is indeed well known that gauge modes of the Maxwell 
field obey a coupled set of eigenvalue equations. In arbitrary gauges, one 
cannot decouple these modes. This can be proved by trying to put in diagonal 
form the 2x2 operator matrix acting on the modes b\ and q. In our problem, 
however, with our choice of gauge-averaging functional and gauge parameter, 
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gauge modes can be disentangled, and a simpler method to achieve this exists. 
For this purpose, we point out that, since the background is flat, if gauge 
modes can be decoupled, they can only reduce to linear combinations of 



MO = — /p— > ( 31 ) 



Bessel functions, i.e., 



and 

q,(r) = C{v)B v {Mt)^. (32) 

With our notation, C(y) is some constant depending on u, which is obtained 
in turn from I. To find v and C(u), we insert the ansatz ([31]) and ([32]) into 
the system of equations (p9|) and (|30|) , and we require that the resulting 
equations should be of Bessel type for B u (Mr), i.e., 



d 2 Id ,,2 A ^ 
+ - — — M 1 + 



B,(Mr) = 0. (33) 



cZr 2 r dr \ M 2 r 2 

This leads to two algebraic equations for v 2 . By comparison, one thus finds 
an algebraic equation of second degree for C: 

1(1 + 1)C 2 - C - 1 = 0, (34) 

whose roots are C + = j, and C_ = — T^jy- The corresponding values of z/ are 
z/ + = Z - i, and i/_ = Z + |. Hence one finds the basis functions for normal 
and longitudinal perturbations in the interior problem in the form 

I l+3/2 (Mr) Ii-i /2 (Mr) 
k{r) = a h i '—= h a 2 ,i , (35) 



c '(0 = -pjJw/2W^+ ^7 M/2 (Mr)^, (36) 



whereas, from Eq. (EH), transverse modes read 

Ti(r) = / i+ i /2 (Mr)^. (37) 



Last, but not least, ghost modes obey an eigenvalue equation analogous to 
'" 2 + 2 -±- Hi±>l] e, = M%, (38) 



dr 2 r dr r 2 



-/(/•) = ^ i+ i/ 2 (Mr). (39) 



and hence they read 



In the exterior problem, i.e., for r greater than the two-sphere radius R, one 
has simply to replace the modified Bessel functions of first kind in Eqs. (]25f) , 
fl35|)-(|37l) and (|39| ) by modified Bessel functions of second kind, to ensure 
regularity at infinity. 

In our problem, ghost modes are of course decoupled from the modes 
for the electromagnetic potential which occur in the expansions (p0|)-([2^). 
Nevertheless, this does not mean that they do not play a role in the Casimir 
energy calculation. By contrast, we already know (see comments after (|10D ) 
that boundary conditions on the ghost are strictly necessary to ensure gauge 
invariance of the boundary conditions on the potential. It is then clear that 



such boundary conditions, combined with the differential equation (|38D , lead 
to a ghost spectrum whose contribution to the Casimir energy can only be 
obtained after a detailed calculation (e.g. Green-function approach, or 
function regularization). It should not be surprising that ghost terms are 
important, since one has already included effects of other degrees of freedom 
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which should be compensated for M. This issue is further clarified by the 
analysis of the original Casimir problem: two perfectly conducting parallel 
plates. In a covariant formalism, one has to consider the energy-momentum 
tensor for ghosts, which is found to give a non-vanishing contribution to 
the average energy density in vacuum. After taking into account boundary 
conditions for the potential entirely analogous to our Eqs. and (|T8l) , 

one then finds a renormalized value of the zero-point energy in complete 
agreement with the result first found by Casimir. 

By virtue of Eq. (0), the modes a;(r) obey homogeneous Dirichlet condi- 
tions: 

KOW = 0,VZ>0. (40) 
Moreover, Eq. ( |i~9|) implies that the modes bi obey the boundary conditions 

= 0,V/>0. (41) 

dM 

Last, the modes q and TJ, being the tangential modes, obey Dirichlet bound- 
ary conditions (cf. Eqs. (||) and (|9])) 

[d(r)W = [r,(r)] fiM = 0,VZ>l. (42) 

On taking into account how the modes are expressed in terms of Bessel 
functions (see Eqs. (|31)|)-(|3"T|)), one thus finds five sets of eigenvalue 

conditions for the interior and exterior problems, respectively: 

(i) Temporal modes: 

I t+1/2 (MR) = 0,Vl>0, (43) 
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K l+l/2 {MR) =0,VZ >0. 



(44) 



(ii) Decoupled normal mode: 

~ d_ 
dr 

d_ 

dr 



r 3 / 2 K 3/2 (Mr) 



r=R 



o, 



= 0. 



Jr=R 



(45) 
(46) 



(iii) Coupled longitudinal and normal modes (here v = I + 3/2): 

(i/ - 1/2)I' U (MR)I^ 2 (MR) + 3(1/ - l)M^/,- 2 (Mi2) 



Mi? 



+ (i/ - 3/2)I'„_ 2 (MR)I v (MR) = 0, 



(47) 



(z/ - l/2)K' v {MR)K v _ 2 {MR) + 3(i/ - l) Ku ^ R) K U _ 2 (MR) 



MR 



+ (u - Z/2)K' V _ 2 {MR)K V {MR) = 0. 



(48) 



(iv) Transverse modes: 



7, +1/2 (Mi2) = 0,VZ>l, 
K l+1/2 (MR)=0,Vl > 1. 



(49) 
(50) 



(v) Ghost modes (multiplying their ^-function by -2): 



J m/2 (Mi?) = 0,V/>0, 
K l+1/2 {MR) =0,VZ > 0. 



(51) 
(52) 
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The eigenvalue conditions ( 4"5|) and ( f46| ) can be re-expressed in the form 

h /2 (MR) = 0, (53) 
K 1/2 (MR) = 0. (54) 



It is thus clear that, by construction, the contribution of Eqs. (|53D and (Q) 
to the Casimir energy of a conducting spherical shell cancels exactly the 
joint effect of Eqs. ©, flUD, (HD> ©> (H and (H), bearing in mind the 
fermionic nature of ghost fields. In general, each set of boundary conditions 
involving a set of positive eigenvalues {A^} contributes to the Casimir energy 
in a way which is clarified by the (^-function method, because the regularized 
ground-state energy is defined by the equation (for Re(s) > s = 2) 



E, 



which is later analytically continued to the value s = in the complex-s 
plane. Here fi is the usual mass parameter and ( is the (^-function of the 
positive-definite elliptic operator B with discrete spectrum {\k}'- 

( B (s) = Tr L2 (B~ s )=^k)- s - (56) 

k 

In other words, the regularized ground-state energy is equal to — |Cb(— 1/2), 
at least for the cases where Cb( s ) has no pole at s = —1/2, as is the case 
for the problem considered here. In general, however, this will not be true 
but instead one has Res Cs( — 1/2) ~ b 2 , and as a result an ambiguity for 
the ground-state energy proportional to the heat-kernel coefficient b 2 of the 
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operator B remains. In these renormalization procedure has to be 

performed to eliminate this ambiguity. 

Although the eigenvalues are known only implicitly, the form of the func- 
tion occurring in the mode-by-mode expression of the boundary conditions 
leads eventually to Eq(0). The non-trivial part of the analysis is represented 
by Eqs. (|47D and (0). These are obtained by imposing the Robin bound- 
ary conditions for normal modes, and Dirichlet conditions for longitudinal 
modes. To find non-trivial solutions of the resulting linear systems of two 
equations in the unknowns aij and a 2) i, the determinants of the matrices of 
coefficients should vanish. This leads to Eqs. ([47]) and (|8|). At this stage, it 
is more convenient to re-express such equations in terms of Bessel functions 
of order I + 1/2. On using the standard recurrence relations among Bessel 
functions and their first derivatives, one thus finds the following equivalent 
forms of eigenvalue conditions: 



I l+1/2 (MR) 



I' l+l/2 (MR) + ^I l+1/2 (MR) 



K l+l/2 {MR) K[ +l/2 {MR) + -±—K l+1/2 (MR) 



= 0,VZ>1, (57) 

0,V/>1. (58) 



H ^ 1V±1L) ^2MR J 
Thus, the contribution of the coupled normal and longitudinal modes splits 

into the sum of contributions of two scalar fields obeying Dirichlet and Robin 
boundary conditions, respectively, with the I = mode omitted. This cor- 
responds exactly to the contributions of TE and TM modes (Eqs. (|3J) and 
(§)), and gives the same contribution as the one found by Boyer 0]. 

Following our initial remarks, it is now quite important to understand 
the key features of the Casimir-energy calculations in other gauges. For this 
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purpose, we consider a gauge-averaging functional of the axial type, i.e. 

= (59) 

where n M is the unit normal vector field n M = (0,1,0,0). The resulting 
gauge-field operator is found to be, in our flat background, 

P v» = g + V ^ V u + -n'W. (60) 

Note that, unlike the case of Lorenz gauge, the a parameter is dimensionful 
and has dimension [length] 2 . Now we impose again the boundary condition 
according to which the gauge- averaging functional should vanish at DM: 

MA)] 9M = [n^A^ dM = [A r ) dM = 0, (61) 

which implies that all bi modes vanish at the boundary (cf. Eq. ([H])). 

A further consequence of the axial gauge may be derived by acting on the 
field equations 

P^Ay = (62) 

with the operations of covariant differentiation and contraction with the unit 
normal, i.e. 

V^ V A V = 0, (63) 
n^ v A v = 0. (64) 
Equation (|63|) leads, in flat space, to the differential equation 

^ + -A r = 0. (65) 
or r 
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This first-order equation leads to bi modes having the form 



bi = % (66) 



Thus, by virtue of the boundary condition (|B"TD , the modes bi vanish every- 
where, and hence A r vanishes identically if the axial gauge-averaging func- 
tional is chosen with such boundary conditions. 
Moreover, Eq. (|4]) leads to the equation 

■ dAt , 1 9 a \ k n <a^\ 
~ iUJ -fr + ^fr A k =0 ' ( 6? ) 

which implies, upon making the analytic continuation uj — > iM, 

, r dai 1(1 + 1) da „ ,„„ s 
M— - - y - ' — - = 0. 68 
ar r z dr 

At this stage, the transverse modes obey again the eigenvalue equation 
(|28|), whereas the remaining set of modes obey differential equations which, 
from Eq. fl6"2|), are found to be 



d 2 + 2d_ 1(1 + 1) 



dr 2 r dr r 2 



ai - M^±Hq = 0, (69) 



^ _ Mat - M 2 q = 0, (70) 

a = 0. (71) 

In particular, Eq. ( [TT]) is obtained from Eq. (|70f ) (when 1 = 0), which is 
a reduced form of the Eq. PkyA v = upon bearing in mind that bi modes 
vanish everywhere. 
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Last, but not least, one should consider the ghost operator, which, in the 
axial gauge, is found to be 

This leads to ghost modes having the form 

si = e ,i e- Mr . (73) 



On the other hand, following the method described after Eq. (|T0|), one 
can prove that, also in the axial gauge, the ghost field should vanish at the 
boundary, to ensure gauge invariance of the whole set of boundary conditions 
on the potential. It is then clear, from Eq. (|73|) , that ghost modes vanish 
everywhere in the axial gauge. 

On studying the system (|5^)-([7D|) one has first to prove that these three 
equations are compatible. This is indeed the case, because differentiation 
with respect to r of Eq. (|68D leads to a second-order equation which, upon 
expressing from Eq. Q6"8| ) and from Eq. (|7(J) , is found to coincide with 
Eq. (69). Thus, we have a system of two second-order differential equations 



for two functions ai and q. However, these functions are not independent, 
in that they are connected by Eq. (0). Hence for every value of / one 
has one degree of freedom instead of two. Finally, we have for every I two 
degrees of freedom, one resulting from Eqs. (|6^)-([70|), and another, i.e. the 
transverse mode T). Thus, an estimate of the number of degrees of freedom 
which contribute to the Casimir energy coincides with that in other gauges. 
Moreover, the parameter a does not affect the Casimir energy, since a does 
not occur in any of the eigenvalue equations. 

16 



Unfortunately, we cannot obtain the exact form of the solutions of Eqs. 
( |69|) and (|70|) in terms of special functions (e.g. Bessel or hypergeometric) . 
This crucial point can be made precise by remarking that, if it were possi- 
ble to disentangle the system (|69|) and (ffOl), one could find some functions 
ai, /3i, Vi, Wi such that the 2x2 matrix 

/ 1 V l \(A l B,U1 a,\ 

\Wi i)\Q A/VA i J 

has no non- vanishing off-diagonal elements, where the operators Ai,Bi,Ci, L>i 
are the ones occurring in Eqs. (|69| ) and (|70l). For example, the first off- 
diagonal element of such matrix is the operator Aiai + Bi + Vi(Ciai + L>i). 
On setting to zero the coefficients of Jtp and 4- one finds that on = —Vi = 
— , where «cm ^ s a constant. But it is then impossible to set to zero the 
"potential" term of this operator, i.e. its purely multiplicative part. 

Nevertheless, there is some evidence that the axial gauge may be con- 
sistently used to evaluate the Casimir energy. For this purpose we find it 
helpful to consider a simpler problem, i.e. the Casimir energy in the axial 
gauge for the case of flat boundary. In this case the basis functions are plane 
waves 

Ai = A 0)IM e i ( k * x+k * v+k ' z - ut) , (74) 

where the admissible values of k x ,k y ,k z are determined by the boundary 
conditions, which are here taken to be analogous to the case of a curved 
boundary. Let us choose the conducting boundaries parallel to the x- and 
y-axes, while the vector n M is directed along the z-axis. Then the condition 
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of compatibility of field equations in the axial gauge is reduced to 

/ —k 2 ujk x loky ujk z 

jj wh x ky-\-k 2 ui 2 k x ky k x k z ^ 

(75) 

Direct calculation shows that this determinant is equal to 



/ -k 2 


Ldk x 


ujky 


tok z 




uk x 


k 2 + k 2 z -co 2 


\\jy 


k x k z 


uiky 


k x ky 


kl + k 2 - u 2 


kyk z 




\uk z 


k x k z 


hy h z 


k 2 x + k 2 -u 2 ^ 





h 2 z (u 2 -k 2 )\ (76) 



D = - 

a 

Hence we have reproduced the correct dispersion relation between energy u 
and wave number k (to every admissible value of k there correspond two 
contributions to the Casimir energy of the form uo =\ k \). Of course, no 
non- vanishing ghost modes exist, once that the axial type gauge-averaging 
functional is set to zero at the boundary (cf. (161]), (|7"i2| ) and (|T3])). Inter- 
estingly, on imposing the boundary conditions, one can set to zero the A z 
component of the electromagnetic potential as was done with A r in the spher- 
ical case. In this case the compatibility condition of field equations is reduced 
to the vanishing of the determinant of a 3 x 3 matrix, obtained by omitting 
the fourth row and the fourth column in Eq. ([75]). One can easily see that 
the determinant of this 3x3 matrix coincides with that of the 4x4 matrix 
up to a multiplicative factor -, which does not affect the dispersion relation. 

To sum up, a complete correspondence can be established between the 
key features of the axial gauge in the cases of flat and curved boundary: 
the a parameter does not affect the Casimir energy, the component of the 
potential orthogonal to the boundary vanishes, ghost modes vanish, and 
only two independent degrees of freedom contribute. On going from the 



flat to the curved case, however, the analysis of the dispersion relation is 
replaced by the problem of finding explicit solutions of Eqs. (|68|) -(fr0|), with 
the corresponding eigenvalue conditions. This last technical problem goes 
beyond the present capabilities of the authors, but the exact results and 
complete correspondences established so far seem to add evidence in favour 
of a complete solution being in sight in the near future. 

We have studied an approach to the evaluation of the zero-point energy 
of a conducting spherical shell which relies on a careful investigation of the 
potential and of ghost fields, with the corresponding set of boundary condi- 
tions on An perturbations and ghost modes. When Boyer first developed his 
calculation, the formalism of ghost fields for the quantization of gauge fields 
and gravitation had just been developed, and hence it is quite natural that, 
in the first series of papers on Casimir energies, the ghost contribution was 
not considered, since the emphasis was always put on TE and TM modes 
for the electromagnetic field. On the other hand, the Casimir energy is a 
peculiar property of the quantum theory, and an approach via path-integral 
quantization regards the potential and the ghost as more fundamental. This 
is indeed necessary to take into account that Maxwell theory is a gauge the- 
ory. Some of these issues had been studied in the literature [§, |J, including 
calculations with ghosts in covariant gauges, but, to our knowledge, an ex- 
plicit mode-by-mode analysis of the and ghost contributions in problems 
with spherical symmetry was still lacking in the literature. The contributions 
of our investigation are as follows ||. 
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First, the basis functions for the electromagnetic potential have been 
found explicitly when the Lorenz gauge-averaging functional is used. The 
temporal, normal, longitudinal and transverse modes have been shown to 
be linear combinations of spherical Bessel functions. Second, it has been 
proved that transverse modes of the potential are, by themselves, unable to 
reproduce the correct value for the Casimir energy of a conducting spheri- 
cal shell. Third, it is exactly the effect of coupled longitudinal and normal 
modes of which is responsible for the value of AE = 0.09hc/2R found 
by Boyer. This adds evidence in favour of physical degrees of freedom for 
gauge theories being a concept crucially depending on the particular problem 
under consideration and on the boundary conditions. Fourth, ghost modes 
play a non-trivial role as well, in that they cancel the contribution resulting 
from transverse, decoupled and temporal modes of the potential. Fifth, the 
axial gauge-averaging functional has been used to study the Casimir energy 
for Boyer's problem. Unlike the case of the Lorenz gauge, ghost modes and 
normal modes are found to vanish, and one easily proves that the result is 
independent of the a parameter. A complete comparison with the case of flat 
boundary has also been performed, getting insight into the problems of inde- 
pendent degrees of freedom and non-vanishing contributions to the Casimir 
energy. 

Indeed, recent investigations of Euclidean Maxwell theory in quantum 
cosmological backgrounds had already shown that longitudinal, normal and 
ghost modes are all essential to obtain the value of the conformal anomaly 
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and of the one- loop effective action Further evidence of the non-trivial 
role played by ghost modes in curved backgrounds had been obtained, much 
earlier, in Ref. [§]. Hence, we find it non-trivial that the ghost formalism 
gives results in complete agreement with Boyer's investigation. Note also 
that, in the case of perfectly conducting parallel plates, the ghost contribution 
cancels the one due to tangential components of the potential (see Sec. 4.5 
of Ref. 0]). This differs from the cancellations found in our paper in the 
presence of spherical symmetry. 

The main open problem is now the explicit proof that the Casimir energy 
is independent of the choice of gauge-averaging functional. This task is as 
difficult as crucial to obtain a thorough understanding of the quantized elec- 
tromagnetic field in the presence of bounding surfaces. In general, one has 
then to study entangled eigenvalue equations for temporal, normal and longi- 
tudinal modes. The general solution is not (obviously) expressed in terms of 
well known special functions. A satisfactory understanding of this problem 
is still lacking, while its solution would be of great relevance both for the 
foundations and for the applications of quantum field theory. In covariant 
gauges, coupled eigenvalue equations with arbitrary gauge parameter also 
lead to severe technical problems, which are described in Ref. ||. A fur- 
ther set of non-trivial applications lies in the consideration of more involved 
geometries, where spherical symmetry no longer holds, in the investigation 
of media which are not perfect conductors, and in the analysis of radiative 
corrections 0, [H]]. Thus, there exists increasing evidence that the study of 
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Casimir energies will continue to play an important role in quantum field 
theory in the years to come, and that a formalism relying on potentials and 
ghost fields is, indeed, crucial on the way towards a better understanding of 
quantized fields. 
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